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Abstract 

Within the frameworks of the hght-front quark model (LFQM) and chiral perturbation theory 
(ChPT) of 0{p^), we reevaluate the form factors of the — > 7 transition. We use these form 
factors to study the decay of e'^u^'j, which is dominated by the structure dependent con- 

tribution. We show the differential decay branching ratio as a function of x = 2E^/mK, where 
(mx) is the photon energy (kaon mass). Explicitly, we find that, in the standard model with 
the cut of X = 0.01 (0.1), the decay branching ratio of — > e^Uej is 1.54 (1.44) x 10^^ and 
1.57 (1.47) X 10-^ in the LFQM and ChPT, respectively. 
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I. INTRODUCTION 



Experimentally, both decays of e^Ve and /i^z/^ have been precisely measured with 

the decay branching ratios being (1.55±0.05) x 10"^ and (63.44±0.14) x 10~^ jl|, respectively. 
The smallness of the electron mode can be easily understood as it is helicity suppressed with 
the suppression factor of m^/m^ ~ 2 x 10~^ in comparison with the muon mode. For the 
corresponding radiative decays of —>■ i^i'i'y {i = e,fi), it is known that they receive two 
types of contributions: "inner bremsstrahlung" (IB) and "structure-dependent" (SD) 
For the decay of K'^ —>■ e^Veli while the IB contribution is still helicity suppressed and 
contains the electromagnetic coupling constant a as well, the SD part gives the dominant 
contribution to the decay rate as it is free of the helicity suppression. Similarly, the SD 
contribution is also important to the decay of — > yU~^z/^7 \a\. 

In the standard model (SM), the decay amplitude of the SD part involves vector and 
axial- vector hadronic currents, which can be parametrized in terms of the vector form factor 
Fy and axial- vector form factor F4, respectively. However, the experimental determinations 



on these form factors are poorly given and model-dependent 
experimental results on the decay rate of e^Vel in Ref. 



7[. In particular, the 
, 0, 0] were based on 



the assumption of Fy and Fa being some constant values in the chiral perturbation theory 
(ChPT) at O(p^) ^. In the ongoing data analysis of the E949 experiment at BNL, more 
precision measurements on the decay of K'^ e^i^el are expected ^ and thus, the model- 
independent extractions of the SD form factors are possible. Theoretical calculations of Fy 
and -p4_in the K'^ 7 transition have been previously done in the ChPT at 0{p^) ^ and 



0{p ) [9|, [lOj. However, the results of the ChPT at 0{p ) [10] have not been fully applied to 



the decay of — > e"'"z/e7 yet. Moreover, it is important if we could obtain information on 
Fy^A in some QCD model other than the ChPT. For this purpose, in the present study we 
will also evaluate Fy^A in the light front quark model (LFQM) . We will use the form factors 
in both ChPT and LFQM to examine the decay of —>■ e'^i'e'j. 

This paper is organized as follows: We present the relevant formulas for the matrix 
elements and form factors for the decay of K'^ e^i'el in Sec. II. In particular, we study 
the transition form factors of — 7 in the ChPT of 0{p^) and LFQM. In Sec. HI, we 
describe the differential decay rate of e+z/e7- In Sec. IV, we show our numerical 

results on the form factors and the decay branching ratio in both ChPT and LFQM. We will 
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also illustrate the differential decay branching ratio as a function of x = 2E^/mK, where 
and ttik are the photon energy and kaon mass, respectively. We give our conclusions in 
Sec. V. 



II. MATRIX ELEMENTS AND FORM FACTORS 



In the SM, the amplitude of the decay —>■ e^Vel (-^627) can be written in terms of IB 



and SD contributions, given by 
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M = MjB + MsD, 

MjB = ie—=sin9cFKmee*^K°', 
v2 

MsD = -ie'^sine,elL,m\ (1) 
v2 



where 



K = U{p^){l + 75 V{pe), 

\Pk -q 2pe ■ g / 

= u(p^)7^(l -75)t;(pe), 

H^-" = — i-g'^PK ■ q + pW) + ^—e'"'''^q.PKp , (2) 
rriK rriK 

ea is the photon polarization vector, px, Pu, Pe, and q are the four- momenta of K~^, u^, e"*", 
and 7, and and Fy^(v) are the K meson decay constant and the axial-vector (vector) 
form factor corresponding to the axial-vector (vector) part of the weak currents, defined by 

{0\sY15u\K+{pk)) = -iFkP^, 
{^{q)\u^^^,s\K{pK)) = [{p ■ q)e*^ - {e* ■ p)q^] , 

{^iq)\uYs\KipK)) = ze^e'^^^'^e.^g^p,, (3) 

"^ir- 
respectively, with p = Pk — q being the transfer momentum. We note that AiiB in Eq- ^ 

is suppressed due to the small electron mass rrie. In the decay of — > e^u^'y, the form 

factors FAy in Eq. are the analytic functions of p"^ = {pk — qY in the physical allowed 

region, given by 

ml<p^ <m\. (4) 

In the following discussion, we will first summarize the formulas for Fy^A in the ChPT 
and then study these form factors in the LFQM. We note that similar calculations for the 
P — » 7 (P = K^, D, B) transitions in the LFQM have been performed in Refs. 



A. Chiral Perturbation Theory 



The chiral Lagrangians contain both normal and anomalous parts. At orders of p™, the 



non-anomalous and anomalous Lagrangians of and £a™'' relevant to the K'^2'y decay are 
given by 



,{m) 



(2) 



Li [Tr{D^UD^'U^)Y + L2Tr{D^UDM^)Tr(D^'UD''U^) 
+L-iTr{D^UD^'U^D^UD''U^) + UTr{D^U D^'U^)Tr{xU^ + Ux^) 
+L,Tr{D^UD^U\xU^ + Ux^)) + L, [Tr{xU^ + f/x^)]' 
+L7 [Trix^U - U^x)f + L,Tr{xU\U^ + Ux^Ux^) 



4'^ = yi7(X+/iM./^'^^)+Z/l8(X+)(V/i''^)+l/8l(X+/+M^/r) 

+iy%-,{f\f,vu^'x+u'') + iyim{f-t,v[f-\ Zip]) + 1/102 (x+Z-M^/-'') 

+yU)^{X+){f-iJ.uf-) + 2/104 X-]) + yi09(Vp/-M^ v'' 

+2/ll0(Vp/+M.[/^'^'',w1) + --, 



and 



167r2 



167r2 

4^) = iC,e^''"^{x-U,uU.p)+iCr,e>^''''^{x-[f+,uJ-.^]) 
+ C22e^^"^(w~{V7/+7^>/+"/3}) + --, 



(5) 



(6) 



respectively, where F is the meson decay constant in the chiral limit, Lj, yj and Ck are 
unrenormalized coupling constants, U is the unitary matrix, parametreized by 

.V2 



U = exp 



v 



TT 

K- 



7T^ , JL 1^0 



(7) 
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Lf^iy and are the field-strength tensors of external sources, given by 



L 
R 



(8) 



and the definitions of all other fields can be found in Ref. 



io| 



From the chiral Lagrangians in Eqs. ([5]) anc 
tions to Fv at O(p^) for the K'^2'y decay to be 

256 



6l) , one obtains the tree and loop contribu- 



io| 



rriK 
AV2tt^Fk 
1 



1 



64 



Iml In f In f ^ U 3mi In ' ^ 



2 / [xm^ + (1 - a;)mj^ - x(l - x)p J In ( ^ — ] dx 



2 / [xml + {l 



xjm 



K 



.xml + {l-x)m\-x{l-x)p^' 
x\l — x)p J In I ■ I ax 



— 4 / mt In 



(9) 



where the wave function and decay constant renormalizations have been included and 
are the renormalized coefficients. From Eq. ([5]), the tree and loop contributions to F4 of 
0{p^) lead to |lo| 



Fk 
rriK 



rriK 



6F|(27r)^ 



:[U2.65{mj^ - p^) - 198.3] 



m: 



{4L; + 7LI + lL[,)mi In -| + 3 (L^ + L\,)m' In 



/ mt 



4V2Fj, 

+2 {8L[ - 4LI + ALl + 71^ + lL\^)m], In 
4v/2mi^ 



3F| 



{2ml{lSy\, - 2yl, - 6y^2 + 2yl, + 3y^4 - yl, + 6y 



103/1 



+2m^(18?/[7 + 36yl8 - %8i - 12i/^2 + 4l/83 + 6z/84 + %85 - 3y 



100 



(10) 



where L[ and ?/[ are the renormalized coupling constants. Note that the first terms in Eqs. 



([9]) and ffTOj) correspond to Fy and F4 at 0(p'^) [4|, |l8|, respectively. We remark that the 
expressions of Eqs. ([9]) and (fTOl) have not been explicitly shown in the literature [sl, lol |. 
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B. Light Front Quark Model 



In the framework of the LFQM |13l . Il4j . [iSj, the physical accessible kinematics region 
is < < M]^ due to the time-like momentum transfers. The general structure of the 
phenomenological light front (LF) meson wave function is based only on the qq Fock space 
sector. It can be expressed by an anti-quark s and a quark u with the total momentum 
(p + q) such as: 

\Kip + q)) = J2 [[dh][dk2]2{2TTf5\p + q-h-k2) 

X ^'^^^-{z,k^)bt{k,,X,)d+{k2,X2m, (11) 



A1A2 



where ^^'^^ is the amplitude of the corresponding s{u) and fci(2) is the on-mass shell LF 
momentum of the internal quark. The LF relative momentum variables {z, k±) are defined 
by 

k^ = zi{p + qy, k:^ = Z2ip + qy, Zi + Z2 = l, 

ki± = Zi{p + q)± + k±, k2± = Z2{p + q)±- k±, (12) 

and 

<i>'^''{z,k^)=( ^j^^ ^\ u{k,,\,)^\{k2A2)<P{z,k^), (13) 

with 0(2;, fcx) being the space part of the wave function, which depends on the dynamics. 
The amplitude of 0(2;, k^ can be solved in principles by the LF QCD bound state equation 



19l . l20[ |. However, we use the Gaussian type wave function in this study: 



K 



(,.,j.A,^/_^exp -^1. (14) 



From Eqs. f|TT|) - f|T4|) . the hadronic matrix elements in Eq. ([3]) are found to be 
{7(9) 157" (1 - 75) «l K(p + 1))= j ^A/t 



' 7 '2 9 I • J 2 9 I • / V 10) I '2 9 I • 

k2 — + le kf — m„ -|- te k{' — + le 

+ {u^s,k[{k^)^k'2{k2))\, (15) 
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where Kk is a vertex function related to the us bound state of the K meson, k2 = q — ki 
and k'^ = [p + q) — k2 = ki + p. After integrating over the LF momentum k^ in Eq. f[T5l) . 
we get 

{^iq)\sYil-l5)u\K{p + q)) 

= r\d%]\ \_ +i^--s,k[{k,)^k',{k,))] , (16) 



where 



2{2T:y^k{^k^kt 

/^Hfc- = Tr{ 75(-^2 + mu)ieui{Hx + mu)l^{\ - 75)(^'i + m,) 

Ion ' 



2 ^ /j^? 

^ion — 5 ^1(2) ~ Von ~ ^2(1)011 ^ ^1 — Q ~ ^2ori 5 (I'i^) 

with {on} representing the on-shell particles. For the kaon, the vertex function Ap in Eqs. 



ffTSjl and ([TID is given by |2l|, l22|: 



^^'^^'"^ ' k ) (18) 

To calculate the right hand part of Eq. f|T6l) . we choose a frame with the transverse momen- 
tum = so that p^ = p^p^ > covers the entire range of the momentum transfers. Here, 
we have used the LF momentum variables {x,k±). Hence, the relevant quark momentum 
variables in Fig. 1 are 

k'+ = (1 - z'){p + g) + , = z'{p + g)+, k[^ = (1 - z')q^_ + k'^, = z'q^_ - /c^ , 

A;+ = (1 - z)q+, k^ = zq+, ki± = (1 - z)q± + k±, k2± = zq± - k± . (19) 

By considering the good component as "/x = +", the hadronic matrix elements in Eq. ([3]) 
can be rewritten as: 

{l{q)\sXl5U+\K{p + q)) = -e-^(el ■gx)p+, 

ZITLk 

(7(g)MK|K(p + g)) = -te^e'^eU.p-'. (20) 



Using Eq, (1191) . the trace part I^'^ in Eq. f[T7j) can be carried out. By comparing Eq. ([16 
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with Eq. (l20l) . we obtain the form factors Fv,a to be: 



where 



3 ml + k'j_ 3 ml + k\ 

2 — 2;' (m^ — fc^B 1 + (1 — 2;')(ms — m„)A;^B 

3 ml + A;2 3 + ^ 1 



A = (1 — 2^)2;' (ms — m„) — 2zmu , 

i? = (1 — 2z)2;'ms + + (1 — 22;)(1 — z )mu , 
, ,2, / -2(1 [dk. 



e 



2[M2 - (m, - m„)2]y V rf^ V 2a; 
1 d^{z, k 



^{z,k\) dk 



2 ' 



p2 



(21) 



^' = ^ (^1 - ^ J , k = {k^, k,) , 

. / \ 4 / 1\ , , ml - ml 

Af„^ = M^ + iI±^. (22) 



III. DIFFERENTIAL DECAY RATE 

In the rest frame, the partial decay rate for e^Vel is given by 

dV = 777^^ I M p dE,dE, , (23) 

where ii^^ and E'e are photon and electron energies, respectively. To describe the kinematics 
of e^Veli we introduce two dimensionless variables, defined by x = 2E^/mK and 

y = 2Ee/mK, with their physical allowed regions being 

< X < 1 — Ve , 
1-X + < y < 1 + re , (24) 

1 — X 
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where Ve = m\/m\;. The relation between the transfer momentum and x is given by: 

/ = m^(l-x). (25) 



From Eqs . ([T]) and fl23l) . we obtain the double differential decay rate of e^Vel as 



m 



dx dy 647r^ 
where A = {x + y — 1 — re)/x and 



(26) 



A{x,y) = AiB{x,y) + AsD+{x,y) + AsD-{x,y) + AiNT+{x,y) + AiNT-{x,y) 



AiB{x,y) 



K\ 



mj^Xx'^ 



x^ + 2(l - re) 1 - a; 



A 



, ,1 — X — - 

1 - A V A 

AsD-i^,y) = \Fv-FA\\^iy-\), 

At / 

AiNT+{x, y) = -Re[FK{Fv + Fa)*] (l - x - 

rriK V 

AjNT-{x,y) = ^Re[FK{Fv - FaT]^ ~ ^ ^ ^ 



A 



rriK A 
By integrating out the y variable in Eq. fl26|) . we obtain the differential decay rate as a 

function of x to be 



where 



A{x) = AjBix) + AsD+ix) + AsD-ix) + AjNT+ix) + AiNT-ix) 



(28) 



Aib{x) 



AreF^ \{x + re- l)[x'^ + 4(1 - re)(l - x)] 



m 



K 



AsD+ 


(x) 




AsD- 


(x) 


= \Fy-FA?x' 


AlNT+ 


(x) 


= ^Re[FK{Fv 
rriK 


-^INT- 


(x) 


= ^Re[FK{Fv 
rriK 



1 — X 

X^ + 2(1 — re)(l — X + Te 
X 

1 — X Tp 



In 



1 — X 



2 6(l-x)2_ 



1 — X Tp 



+ 



2 6(1 -x) 

1 - X rl , 1 '^e 

-, r + fe In 

2 2(1 -x) 1-x 



-1 + 3x ^ — 2xrp 



+ (x — rS] In ■ 



2 ■ 2(l-x) ' ''""^l-x_ 
It is clear that the contributions to the decay rate from the IB and INT^ parts are sup- 
pressed due to the small electron mass. 



IV. NUMERICAL RESULTS 



The numerical values of FAyip^) in the ChPT of 0{p^) have been shown in Figs. 5 and 



6 of Ref . 



10|. To compare these values with those in the LFQM, we plot the results in Figs. 



[T]and[2J In these figures, we have also included the results in the ChPT at 0{p^). For the 



0.14. 
0.12 

0.08 ■ 
0.06 ■ 
0.04. 
0.02 ■ 
0.00 ■ 



■ChPT 0(p'*) 
■ChPT O(p') 
LFQM 



0.00 0.05 0.10 0.15 

p' (GeV') 



0.20 



0.25 



FIG. 1: Fy(p^) as a function of the transfer momentum 



calculations of the ChPT 



we have taken itik = 0.495, m^r = 0.14, m„ = 0.55 and irip 



0.77, Fk = 0.112 GeV and the renormalized coefficients of (L^, L^, L|, L^, L\q), (Q, C[i, Q2) 
and (i/[oo, Z/[o4, vim. Vlio) to be (0.53, 0.71, -2.72, 6.9, -5.5) x 10-^ [23], (0.013, -6.37, 6.52) x 
l()~^GeV-^ [24] and (1.09, -0.36,0.40-0.52) x 10-7F|. y, respectively. For some other 
possible sets of coefficients, see Ref. 10| as well as the recent review in Ref. 26|. We note 
that we have ignored the contributions from j9^-nondependent terms involving On the 
other hand, the p^-dependence of F4(p^) for the ChPT at 0{p^) are insensitive due to the 



small contributions related to yl [lOj. We emphasize that as illustrated in Figs. [T] and [2l 
the form factors -FV,a at O(p^) in the ChPT are constants [4]. To evaluate the form factors 
of FyA from Eq. fl2Tl) in the LFQM , we have used rriu = 0.26, rris = 0.37 and ujk = 0.382 
in GeV. In Tabled we explicitly display the values of Fv,a{p'^ = 0). 

By integrating out the variable x in Eq. (l28l) . in Table HTl we give the decay branching 
ratio of K+ e+Vel in (a) the ChPT at 0(/), (b) the ChPT of 0{p^) and (c) the LFQM. 
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ChPTOCp") ^\ 

ChPTO(p') '^.^ 
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O.bo ' 005 ' 010 ' 015 ' O20 ' 025 
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FIG. 2: Fa{p'^) as a function of the transfer momentum p^. 

TABLE I: The form factors of Fv(0) and Fa{0) in (a) the ChPT at O(p^) Q], (b) the ChPT of 
0(/) and (c) the LFQM. 



Model 


Fv{0) 


Fa{0) 


(a) 


0.0945 


0.0425 


(b) 


0.082 


0.034 


(c) 


0.106 


0.036 



Here, as the IB term diverges at the limit of x — corresponding to p"^ — > Pmax ~ ^ 
we have used the cuts of x = 0.01 and 0.1, respectively. With the cuts, from Table HTl we 
see that both IB and INT^ contributions are much smaller than the SD^ ones, which are 
insensitive to the cut. We remarks that in Table [TTl our results for the SD"*" contribution 
to the decay branching ratio in the ChPT of 0(»^) and LFQM are 1.15 and 1.12 x 10^^, 
which are smaller than that of 1.52 ± 0.23 x 10"^ fl, quoted by the PDG , respectively. 
Note that the value in the PDG was based on the combination of the data in Refs. [sl and 
6|, in which large constant values of Fa + Fy = 0.1501q;q23 0.147 ± 0.011 were used, 
respectively. It is clear that to compare the data with the theoretical predictions, proper 
form factors should be used in the data analysis. 
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To show the behaves of the various contributions in the ChPT and LFQM, we present the 
IB and SD^ parts of the differential decay branching ratio as functions of x in Fig. [31 Here, 
we do not plot the INT^ contributions in Fig. 4 as they are vanishingly small. As shown in 
the figure, in the small x region there is an enhancement for the IB part, whereas those from 
the SD^ parts are close to zero. In Fig. HI we also display the spectrum of the differential 
decay branching ratio vs. x in the ChPT at both O(p^) and O(p^) and the LFQM. 

TABLE II: The decay branching ratio of e.'^i'el (in units of 10~^) in (a) the ChPT at O(p^), 

(b) the ChPT of 0{p^) and (c) the LFQM with the cuts of x = 0.01 and x = 0.1, respectively. 



Model 


Cut 


IB 


SD+ 


SD- 


INT+ 


INT- 


Total 


(a) 


x = 0.01 
X = 0.1 


1.65 X IQ-i 
0.69 X IQ-i 


1.34 
1.34 


1.93 X 10-^ 
1.93 X IQ-^ 


6.43 X 10-5 
6.43 X 10-^ 


-1.10 X 10-3 
-1.10 X 10-3 


1.70 
1.60 


(b) 


X = 0.01 
X = 0.1 


1.65 X IQ-^ 
0.69 X IQ-i 


1.15 
1.15 


2.58 X 10~i 
2.58 X 10~i 


6.22 X 10-^ 
6.22 X 10-5 


-1.21 X 10-3 
-1.21 X 10-3 


1.57 
1.47 


(c) 


X = 0.01 
X = 0.1 


1.65 X 10^1 
0.69 X 10-1 


1.12 
1.12 


2.59 X IQ-^ 
2.59 X IQ-^ 


4.33 X 10-5 
4.33 X 10-5 


-1.29 X 10-3 
-1.29 X 10-3 


1.54 
1.44 



50-, 



40- 



X 
T3 

^ 30- 

> 

t 

^ 20- 

CO 
T3 



0- 



0.0 0.2 0.4 0.6 0.8 1.0 

x(=2E /m^) 



FIG. 3: The IB and SD^ parts of the differential decay branching ratio as functions of x = 2E^/mK- 
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IB 

SD-HChPT O(p') 

SD- ChPT O(p') 

SD+ LFQM 

SD- LFQM 





FIG. 4: The differential decay branching ratio as a function of x = 2E^/mK- 

From Fig. HJ we see that in the region of x < 0.7 or i?^ < 173 MeV, the decay branching 
ratio in the LFQM is much smaller than that in the ChPT at 0{p^). On the other hand, in the 
region of x > 0.7 the statement is reversed. However, if we only consider the contributions in 
the ChPT at O(p^), the conclusion is weaker. In Table UTTt we illustrate the decay branching 
ratio in the regions of 0.1 < x < 0.7 and 0.7 < x < 1 from the various approaches, 
respectively. The main reasons for the differences are due to the form factors. The form 

TABLE III: The decay branching ratio of — > e^VeJ (in units of 10~^) in the regions of 
0.1 < X < 0.7 and 0.7 < x < 1 from the various approaches, respectively. 



Region 


ChPT of 0{p^) 


ChPT of 0(/) 


LFQM 


O.K X < 0.7 


0.871 


0.871 


0.541 


0.7 < X < 1 


0.733 


0.606 


0.902 



factors of the ChPT at O(p^) are constant and straight lines at 0{p^), whereas in the LFQM 
they are the overlap between the wave functions of the K meson and photon and become 
zero when x ^ or ^ Pmax = ''^x- 1^ is clear in the future data analysis such as the 
one at the experiment BNL-E949 Q , one could concentrate on these two regions to find out 
which model is preferred. 
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V. CONCLUSIONS 



We have studied the axial-vector and vector form factors of the — ^ 7 transition in 
the LFQM and ChPT of 0{p^). Based on these form factors, we have calculated the decay 
branching ratio of e^Vel- We have demonstrated that the SD part gives the dominant 

contribution to the decay in the whole allowed region of the photon energy except the low 
endpoint. Explicitly, we have found that, in the SM with the cut of x = 0.01 (0.1), the 
decay branching ratio of e^Vel is 1.54 (1.44) x 10~^ and 1.57 (1.47) x 10~^ in the 

LFQM and ChPT, respectively. Future precision experimental measurements on the decay 
spectrum js] should give us some useful information to determine the SD contribution as 
well as the vector and axial- vector form factors. 
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